Shear dynamo problem: Quasilinear kinematic theory 
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Large-scale dynamo action due to turbulence in the presence of a linear shear flow is studied. Our treatment 
is quasilinear and kinematic but is non perturbative in the shear strength. We derive the integro-differential 
equation for the evolution of the mean magnetic field, by systematic use of the shearing coordinate transforma- 
tion and the Galilean invariance of the linear shear flow. For non helical turbulence the time evolution of the 
cross-shear components of the mean field do not depend on any other components excepting themselves. This 
is valid for any Galilean-invariant velocity field, independent of its dynamics. Hence the shear-current assisted 
dynamo is essentially absent, although large-scale non helical dynamo action is not ruled out. 

PACS numbers: 47.27.W-, 47.65.Md, 52.30.Cv, 95.30.Qd 



Shear flows and turbulence are ubiquitous in astrophysical 
systems. Recent work suggests that the presence of shear may 
open newpathways to the operation of large-scale dynamos 

ill II mill. We present a theory of dynamo action in a 

shear flow of an incompressible fluid which has random ve- 
locity fluctuations due either to freely decaying turbulence or 
generated through external forcing. Of particular interest is 
the case of non helical large-scale dynamo action in shear 
flows. Several direct simulations show that large-scale fields 
can grow from small seed fields under the combined action 
of non helical turbulence and background shear flow 
However, the interpretation of how such a dynamo works is 
not yet clear One possibility that has attracted much attention 
is the shear-current effect [4], in which extra components of 
the mean electromotive force (EMF) arise due to shear, which 
couple components of the mean magnetic field parallel and 
perpendicular to the shear flow. However there is no conver- 
gence yet on whether the sign of the relevant coupling term is 
such as to obtain a dynamo; some analytic calculations 
and numerical experiments Ij]] find that the sign of the shear- 
current term is unfavorable for dynamo action. Moreover, an- 
alytic calculations treat shear as a small perturbation. We are 
interested here in studying the shear dynamo without such a 
restriction. 

Our theory is 'local' in character: In the lab frame we con- 
sider a background shear flow whose velocity is unidirectional 
(along the X2 axis) and varies linearly in an orthogonal direc- 
tion (the Xi axis). The linear shear flow has a basic sym- 
metry relating to measurements made by a special subset of 
all observers, who may be called comoving observers. This 
symmetry is the invariance of the equations with respect to a 
group of transformations that is a subgroup of the full Galilean 
group. It may be referred to as 'shear-restricted Galilean in- 
variance', or Galilean invariance (GI). We introduce and ex- 
plore the consequences of GI velocity fluctuations; not only 
are these compatible with the underlying symmetry of the 
problem, but they are expected to arise naturally. This has 
profound consequences for dynamo action, because the trans- 



port coefficients that define the mean EMF become spatially 
homogeneous in spite of the shear flow. Systematic use of 
the shearing transformation allows us to develop a theory that 
is non perturbative in the strength of the background shear. 
However, we ignore the complications associated with nonlin- 
ear interactions, hence MHD turbulence and the small-scale 
dynamo; so our theory is quasilinear in nature, equivalent to 
the 'first order smoothing approximation' (FOSA). 

Let (ei, 62, 63) be the unit vectors of a Cartesian coordi- 
nate system in the lab frame, X — (Xi,X2,X3) the posi- 
tion vector, and r the time. The fluid velocity is given by 
{~2AXie2+v), where A is the shear parameter and f (X, t) 
is a randomly fluctuating velocity field which is incompress- 
ible (V- V ~ Q) and has zero mean {{v) = 0). The magnetic 
field has a large-scale (mean field) component B{X, r), and 
a fluctuating field, b, with zero mean ((b) = 0). The evolution 
of the mean field is governed by 
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B + 2ABie2 = Vxe + r]V^B (1) 



where £ — (vxb) is the mean EMF. Our goal is to calculate £ 
in terms of the statistical properties of the fluctuating velocity 
field, which we will do using quasilinear theory. This means 
solving the equation for b by dropping terms that are quadratic 
in the fluctuations. We also drop the resistive term, assuming 
that the correlation times are small compared to the resistive 
timescale. So our theory is applicable when FOSA is valid 
H. Then b obeys 



OT 0X2 



b + 2Abie2 = Vx{vxB) (2) 



It proves convenient to exchange spatial inhomogeneity for 
temporal inhomogeneity, so we get rid of the {Xid / 8X2) 
term through a shearing transformation to new spacetime vari- 
ables, 



Xl 



Xi 



X2=X2+ 2AtXx , X3^X3, t^T (3) 
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We also define new variables, if (a;, = B{X,T),h{x,t) = 
b(X, t) and u{x, t) — v{X^ r), which are component-wise 
equal to the old variables. 
Then equation (|2|i becomes 
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H (4) 



Not only do sheared coordinates get rid of spatial inhomo- 
geneity, but in quasilinear theory the evolution equation (|4|i 
does not contain spatial derivatives of h{x^t). The equations 
for hi and can be integrated directly. The /ii so obtained 
can be substituted in the equation for /12 : there occur double- 
time integrals which can be manipulated to give expressions 
with only single-time integrals, by interchanging the order of 
the integrals. Then the particular solution for h{x, t) is given 
in component form by 



dt'\ 



'•ml 



- 2A{t - t')S,n2u[i] [H[ + 2At'5i2H[] 



dt' [u'l + 2At'Si2u\] [H'„ 



ml 



2A{t - t')6rn2H[i] 

(5) 



where primes denote evaluation at spacetime point {x,t'). 
We have also used notation = (dum/dxi) and Hmi = 
(dH^/dxi). 

The expression in equation ^ for h should be substituted 
in 5 = (vxb) = (uxh). Following standard procedure, we 
allow ( ) to act only on the velocity variables but not the mean 
field; symbolically, it is assumed that (uuH) = (uu) H. 
After interchanging the dummy indices {I, m) in the last term, 
we find that the mean EMF is 

dt' \aa - 2A{t - [H'l + 2At'5i2 H[] - 

I dt' [ 57,„u + 2At'5^m2 mil] [H'l^ - 2A{t - t')6i2 H'^^] (6) 
Jo 

where the transport coefficients {a ,(3 ,r]) are defined in terms 
of the uu velocity correlators by 

ail{x,t,t') = eij,n {uj{x,t)u„ri{x,t')) 
%{x,t,t') = €ij2 {Uj{x,t)uii{x,t')) 

rjirni{x,t,t') = e.ji {uj{x,t)um{x,t')) (7) 

It is physically more transparent to consider velocity statistics 
in terms of the vv velocity correlators, because this is referred 
to the lab frame, instead of the sheared coordinates. By defi- 
nition. 



where X{x, t) = {xi, X2 — 2Atxi, X3) is the inverse of the 
shearing transformation given in equation (|3]l. The velocity 
gradient u^i is 

f d d \ 

Uml = TTTT ~ 2At(5/1 — — Vm = Vml - 2At 6ll Vm2 

\ dXi 8X2 J 

(9) 

where Vmi = {dv,n/ dXi). Using equations ^ and (|9]l in (|7]|, 

ail{x,t,t') = eijm[{vjiX,t)VraliX',t')) 

-2At'Sn {v,{X,t)vm2iX',t'))] 

Puix,t,t') = e,j2[{vj{X,t)vniX',t')) 

-2At'Sn {vj{X,t)vi2{X',t'))] 

%r,u{x,tX) = e,ji{vj{X,t)v^{X' ,t')) (10) 

where the quantities X — {xi , X2 — 2Atxi , 2:3 ) and X' = 
{xi ,X2 - 2At'xi ,2:3). 

We can arrive at some general conclusions for delta- 
correlated-in-time velocity fields. Let the the two-point cor- 
relator taken between spacetime points (il, r) and (il', t') be 
{vi{R,T)vj{R! ,t')) = 5{t ~t')T,j{R,R:,t). We define 
Tiji{R,T) = {dTij / dR'i) The delta-function ensures 

that X and X' occuring in the velocity correlators of equa- 
tion ( [Tol l are equal to each other So (vi{X, t) Vj{X' , i')) — 
5{t - t')T,j{X,X,t), and {v,{X,t)v,i{X' ,t')) = 6{t - 
t') Tiji{X, t). The integrals over time in equation (|6]l can all 
be performed, so the mean EMF is 

— ^ijm\Tjml ^'2At5iiTjm2][Hl + 2AtSl2Hi] 

— £ijl[Tjm + 2AtS,n2Tji] Him (H) 

It is useful to write the EMF in terms of the original variables 
and lab frame coordinates. To this end we transform 



2ATSml TTT^ 1 Bl = Bl,n-2ATSml Bl2 



ax 



dXo 



(12) 

where Bi„i = [dBi/ dXm)- Then the explicit dependence 
of £i on the shear parameter A cancels out, and mean EMF 
assumes the simple form. 



(13) 



Um{x,t) = VmiX{x,t),t) 



(8) 



which is the familiar expression obtained in the absence of 
shear. Thus, shear needs time to manifest and, to see the ef- 
fects of shear explicitly, it is necessary to consider non zero 
correlation times. Henceforth we consider velocity statistics 
with finite correlation times. 

The linear shear flow has a basic symmetry relating to mea- 
surements made by a special subset of all observers. We define 
a comoving observer as one whose velocity with respect to the 
lab frame is equal to the velocity of the background shear flow, 
and whose Cartesian coordinate axes are aligned with those 
of the lab frame. A comoving observer can be labelled by the 
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coordinates, ^ = (^i , ^2, (.3), of her origin at time r = 0. Dif- 
ferent labels identify different comoving observers and vice 
versa. As the labels run over all possible values, they exhaust 
the set of all comoving observers. The origin of the coordinate 
axes of a comoving observer translates with uniform velocity; 
its position with respect to the origin of the lab frame is given 
by 



(14) 



An event with spacetime coordinates {X, t) in the lab frame 
has spacetime coordinates (X, f ) with respect to the comov- 
ing observer, given by 



X = X Xe(T), 



To 



(15) 



where the arbitrary constant tq allows for translation in time 
as well. 



Let 



denote the mean, the 



B(X,f),6(X,f),7}(X,f) 
fluctuating magnetic fields and the fluctuating velocity 
field, respectively, as measured by the comoving ob- 
server. They are all equal to the respective quantities mea- 
sured in the lab frame: f ) , f ) , {^(X , f ) = 

[B{X , t) , b{X, t) , v{X, t)]. That this must be true may 
be understood as follows. Magnetic fields are invariant under 
non-relativistic boosts, so the mean and fluctuating magnetic 
fields must be the same in both frames. To see that the fluctu- 
ating velocity fields must also be the same in both frames, we 
note that the total fluid velocity measured by the comoving 

observer is, by definition, equal to (^—2AXe2 + v{X,f)j. 
This must be equal to the difference between the velocity 
in the lab frame, {-2AXe2 + v{X ,t)), and (-2^^162), 
which is the velocity of the comoving observer with re- 
spect to the lab frame. Using X = X — ^1, we see that 
v{X,f) = v{X,t). Equations ([T]i and are invariant 
under the simultaneous transformations of spacetime coordi- 
nates and fields discussed above. We note that this symmetry 
property is actually invariance under a subset of the full ten- 
parameter Galilean group, parametrized by the five quantities 
(^1, ^2, £,3,To,A); for brevity we refer to this restricted sym- 
metry as Galilean invariance, or simply Gl. There is a fun- 
damental difference between the coordinate transformations 
associated with Gl (equation [TSll and the shearing transfor- 
mation (equation |3]l. The former relates different comoving 
observers, whereas the latter describes a time-dependent dis- 
tortion of the coordinate axes of one observer. Moreover, the 
relationship between old and new variables is homogeneous 
for the Galilean transformation, whereas it is inhomogeneous 
for the shearing transformation. 

Naturally occuring processes lead to G-invariant velocity 
statistics. Let the observer in the lab frame correlate Vi at 
spacetime location (-R, t) with vj at location {R',t'). Now 
consider a comoving observer, the position vector of whose 
origin is given by Xc{t) of equation (fl4l l. An identical exper- 
iment performed by this observer must yield the same results, 
the measurements now made at the spacetime points denoted 



by (R + Xc(t),t) and {R' + Xc{t'),t') in the lab frame 
variables. Therefore, a Gl two-point velocity correlator must 
satisfy the condition, 

{v,{R,r)v,{R\r')) = 

{v,{R + X,(r), r) v,{R! + X,{t'),t')) (16) 

for all {R, R', t, t', |). We also have 

{v,{R,t)v,i{R',t')) = 

{v,iR + X,{t),t) v,i{R' + X,(t'), r')) (17) 

If we now set R = R' = ,T = t ,t' ^ t' and ^ = a;, 
we win have Xc{t) — {xi , X2 — 2Atxi , X3) and Xc{t') = 
{xi ,X2 — 2At'xi ,xs). Therefore Xc{t) and Xc{t') are 
equal to X and X', which are the quantities that enter as ar- 
guments in the velocity correlators of equations ( fTOl ) defining 
the transport coefficients. Hence, (reading equations ( fTSI l and 
( fTTb from right to left), we see that 

{v,{X,t)vj{X',t')) = {v,{0,t)vj{0,t')) ^ R,,{t,t') 
{v,{X,t)vji{X',t')) = {v,iO,t)vji{0,t')) ^ S^ji{t,t') 

(18) 

are independent of space, and are given by the functions 
Rij {t, t') and Siji {t, t'). Symmetry and incompressiblity im- 
ply that i') = Rji{t',t) and Sijj{t,t') = 0. Using 
equations ( fTSI ) in equations ( fTol i. we find that the Gl transport 
coefficients 

(iil(t,t') = eijjn[Sjml{t,t') - 2At' Sii Sjm2{t,t')] 

h [t, t') = ey2 [Sjii (t, t') - 2 At' Sn Sji2 (t, <')] 

%-ml{tX) = e^JlRJ^{t,t') (19) 

are also independent of space. 

Galilean invariance is the fundamental reason that the ve- 
locity correlators, hence the transport coefficients, are inde- 
pendent of space. The derivation given above is purely math- 
ematical, relying on the basic freedom of choice of parameters 
{R, R' , T, t' , ^), but we can also understand the results more 
physically. X and X' can be thought of as the location of 
the origin of a comoving observer at times t and t', respec- 
tively. Gl implies that the velocity correlators measured by 
the comoving observer at her origin at times t and t' must 
be equal to the velocity correlators measured by any comov- 
ing observer at her origin at times t and t'. In particular, this 
must be true for the observer in the lab frame, which explains 
equations ( fTSl l. consequently equations ( fT9l ). We can derive 
an expression for the Gl mean EMF by using equations fT% 
for the transport coefficients in equation (|6]l, and simplifying 
the integrands. Define 



Cjml [t, t ) 
Djm{t,t') 



Sjml{t,t') - 2A{t-t')5m2Sjil(t,t') 
Rjrn{t,t') + 2At'd,n2Rjl{t,t') (20) 
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Then the mean EMF, S{x,t), can be written compactly as 
Jq 

- / dt' [e^Jl - 2A{t - t')Siie,j2] Djm{t,t')H'i„^ 
Jo 

(21) 

The mean field equation ^ for H{x^ t) is 

+ 2^5,2^1 = ( V X + riV^H, (22) 

where (V)^^= d/dXp = {d/dxp + 2Atdpid/dx2). We use 
equation ( I2l1 l to evaluate {'Vx£)^: 

Jo 
Jo 

~2A{t ~ t')S,2 [H[^„, + 2AtS,iH[^„,] } 

(23) 

Equations ( |22] | and ( l23T l form a closed set of integro-diffential 
equations governing the dynamics of the mean field, H{x,t), 
valid for arbitrary values of A. The most visible properties of 
equation ( |23l ) for (V X 5) are: (i) Only the pait of dmi {t, t') 
that is antisymmetric in the indices {i, m) contributes. Indeed 
both Simi and Cimi can vanish for non helical velocity fluctu- 
ations, in which case dynamo action is determined only by the 
Djm terms, (ii) The Dj„i{t, t') terms are such that (Vx5)j 
involves only Hi for i = \ and i = 3, whereas (Vx5)2 de- 
pends on both H2 and Hi . Together with the mean field induc- 
tion equation (|22] | this means that the equations determining 
the time evolution of Hi and Hj, are closed. Thus Hi{x,t) 
(or H^(x^ t)) can be computed by using only the initial data 
Hi{x,Q) (or H^{x,Q)). The equation for H2 involves both 
H2 and Hi, and can then be solved. 

The implications for the original field, JB(X,t), can be 
read off, because it is equal to II{x, t) component-wise (i.e 
Bi{X,T) = Hi{x,t)). Thus, the Djm{t,t') terms do not 
couple either Bi or with any other components, except- 
ing themselves. In demonstrating this, we have not assumed 
that either the shear is small, or that H{x, t) is such a slow 
function of time that it can be pulled out the time integrals in 
equations ( 1211 1 and ( l23T l. Comparing with earlier work (where, 
essentially, both assumptions have been made) we conclude 
that there is no shear-current assisted dynamo of the form 
discussed 0/> where there is explicit coupling of B2 



and Bi in the evolution equation for Bi . Our calculations 
are based on a non perturbative treatment of shear, and this 
makes for a basic departure from earlier work which have 
treated shear perturbatively. Even when the shear is weak, 
two fluid elements which were close together initially would 
be separated by arbitrarily large distances at late times. Thus 
the two-time correlators, which appear naturally in the dy- 
namo problem, have to be handled carefully in the presence 
of shear Moreover, the perturbative treatment of shear is not 
guaranteed to preserve GI, which is a natural and fundamental 
ingredient of our non perturbative approach. 

In conclusion we find that systematic use of the shearing 
coordinate transformation and the Galilean invariance of a lin- 
ear shear flow allows us to develop a quasilinear theory of the 
shear dynamo which, we emphasize, is non perturbative in the 
shear parameter Specifically, we have proved that there is es- 
sentially no shear-current assisted dynamo in the quasilinear 
limit when FOSA is applicable. Moreoever, our results are 
valid for any GI velocity statistics, independent of the forces 
(Coriolis, buoyancy etc) governing the dynamics of the ve- 
locity field. However, large-scale non helical dynamos (i.e. 
with no initially imposed kinetic helicity) are not ruled out, 
and further progress requires developing a dynamical theory 
of velocity correlators in shear flows. 

We acknowledge Nordita for providing a stimulating at- 
mosphere during the program on 'Turbulence and Dynamos'. 
We thank Axel Brandenburg, Karl-Heinz Radler and Matthias 
Rheinhardt for valuable comments. 



* Electronic address: ssridhar@rri.res. in 
^ Electronic address: kandu@iucaa.emet.in' 

[1] A. Brandenburg et al., Astrophys. J., 676, 740 (2008). 

[2] T. A. Yousef et al., Phys. Rev. Lett., 100, 184501 (2008); Astron. 
Nachr. 329, 737 (2008). 

[3] P. J. Kapyla, M. J. Korpi and A. Brandenburg, Astron. Astro- 
phys., 491, 353 (2008); arXiv:0812.1792 (2008); D. W. Hughes 
and M. R. E. Proctor, arXiv:0810.1586 (2008). 

[4] I. Rogachevskii and N. Kleeorin, Phys. Rev. E 68, 036301 
(2003); 70 046310 (2004); Astron. Nachr., 329, 732 (2008). 

[5] A. A. Schekocihin et al., arXiv:08 10.2225 (2008). 

[6] K.H. Radler and R. Stepanov, Phys. Rev E 73, 0563 11 (2006). 

[7] G. Rudiger and L. L. Kitchatinov, Astron. Nachr. 327, 298 
(2006). 

[8] H. K. Moffatt, Magnetic Field Generation in Electrically Con- 
ducting Fluids, Cambridge University Press, Cambridge (1978); 
F. Krause, K.-H. Radler, Mean-field magnetohydrodynamics and 
dynamo theory, Pergamon Press, Oxford (1980); A. Brandenburg 
and K. Subramanian, Phys. Rep. 417, 1 (2005). 



